Metal-superconductor transition at zero temperature: 
A case of unusual scaling 

[ Phys. Rev. Lett. 79, 3042 (1997)] 
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An effective field theory is derived for the normal metal-to-superconductor quantum phase transi- 
tion at T — 0. The critical behavior is determined exactly for all dimensions d > 2. Although the 
critical exponents j3 and v do not exist, the usual scaling relations, properly reinterpreted, still hold. 
A complete scaling description of the transition is given, and the physics underlying the unusual 
critical behavior is discussed. Quenched disorder leads to anomalously strong T c -fluctuations which 
are shown to explain the experimentally observed broadening of the transition in low-T c thin films. 



PACS numbers: 74.20.-z , 64.60.Ak 

There has been much interest lately in quantum phase 
transitions, which occur at zero temperature (T = 0) as a 
function of some nonthermal control parameter. The sub- 
ject of much of this attention have been magnetic tran- 
sitions [Q. A characteristic feature of these transitions 
is that the upper critical dimension dj, above which the 
critical behavior is governed by a simple Gaussian fixed 
point (FP), is lower than that of the corresponding clas- 
sical or finite temperature transition, if any. The reason 
is the coupling of statics and dynamics that is inherent 
to quantum statistical mechanics, which leads to an ef- 
fective increase of the system's dimension from d to d+z, 
with z the dynamical critical exponent. 

In this Letter we study the transition from a normal 
metal to a superconductor at T = 0. We concentrate 
on the realistic case of a superconductor with a nonzero 
density of elastic scatterers, so that the normal state con- 
ductivity is finite. The transition can then be triggered 
by varying either the attractive electron-electron inter- 
action, or the disorder. For such systems we find d+ = 2. 
For d > 2 the transition is governed by a Gaussian FP 
with unusual properties. The resulting critical behavior 
can be summarized as follows. The correlation length £ 
and the order parameter (OP) \& scale like 
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with t the dimensionless distance from the critical point 
(t < in the disordered phase), and the step function. 
The physical energy gap A, as measured in a tunneling 
experiment, is distinct from the OP and scales like 
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The relations (|la|) imply that although the usual critical 
exponents v and (3, defined by £ ~ \t\~" and * ~ \t\P, do 
not exist (y = oo and = oo), one can assign a value to 
their ratio, 



p/v = 2 



(lc) 



in the sense that A £ 2 ~ const. The two-point vertex as 
a function of the wavevector q and the frequency Q, mea- 
sured in suitable units, asymptotically close to criticality 
has the structure 
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Besides the behavior of the correlation length given 
above, one reads off three more critical exponents from 
this result: The exponent r), 7, and z, defined as 
r( 2 )(q,ft = 0) ~ M 2 "'", r( 2 )(q,fi = 0) ~ P, and 
r r ~ £ z , with r r the relaxation time, respectively, are 



77 = 2 



7 = 1 



z = 2 
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where 77 is defined only up to logarithmic corrections Q 
We see that the usual scaling relation 77 = 2 — 7/1/ is 
obeyed (if one puts v = 00 and ignores logarithmic cor- 
rections), while all hyperscaling relations are violated. 

The magnetic penetration depth A and the upper crit- 
ical field H C2 scale like the correlation length and the 
inverse correlation length squared, respectively, 



r 
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The free energy density obeys the scaling law 

/(t,T) = 6- 4 /(iln6,T6 2 ) , (4a) 

with b an arbitrary length rescaling parameter. This 
yields in particular for the specific heat coefficient, 7 = 
d 2 f/dT 2 , the homogeneity law 



7(t,T) = 7 (iln6,T6 2 ) 



(4b) 



This implies, among other things, that the specific heat 
coefficient approaches a constant as the temperature is 



1 



lowered to zero at the critical coupling strength. More 
detailed calculations show a step discontinuity in j(T = 
0) at t = 0. An analogous scaling behavior is found for 
the electrical resistivity, which obeys 



p(t, T) = p(t\nb, Tb 2 ) 



(5) 



Here the step discontinuity leads from a finite value to 
zero. Finally, we find that the fluctuations of the location 
of the critical point are long-range correlated, falling off 
like l/r 2 ' d ~ 2 ). This modifies || the usual Harris criterion 
H for the correlation length exponent to v > l/(d— 2), 
which is obeyed at our Gaussian FP for all d > 2. As 
a result of this long-rangedness, one finds anomalously 
strong fluctuations in the position of the critical point as 
one approaches the quantum FP (although the transition 
remains sharp). For the root-mean square fluctuations 
of the critical temperature we obtain 



AT C /T C ~ Tj d - 2 )/ 2 F(ln(T /T c )) 



(6) 



with To the characteristic temperature of the excitations 
that mediate the superconductivity, and F(x) a func- 
tion that parametrically depends on the disorder. This 
explains the strong broadening of the transition that is 
observed in thin superconducting films at low values of 
T c . Our result is in semi-quantitative agreement with 
recent mesurements of the tunneling conductance G in 
PbBi films [^). In the parameter regime of this experi- 
ment, F(x) ~ x is a reasonable approximation. Defining 
AT C as the width of the region over which G drops from 
80% to 20% of its normal state value, and using d = 2 || 
and T = 300 K, we find that the data of Ref. || agree 
with Eq. (|) to within 20%. 

The unusual scaling properties listed above are in sharp 
contrast to the thermal critical behavior of superconduc- 
tors, which are in the universality class of a d-dimensional 
XY-model. In the remainder of this Letter we derive and 
discuss these results. Our starting point is a very general 
fermionic action that we write in the form 



S = S Q - (K c /itN f ) I dq n c (q) n c {q) 



(7) 



The explicitly written part of this action is the Cooper 
channel interaction term, with K c the appropriate in- 
teraction constant normalized by means of the free elec- 
tron density of states per spin at the Fermi level, Np. 
For an attractive interaction, K c < 0. q = (q, Q n ) 
comprises space and (bosonic) Matsubara frequency la- 
bels, and J dq — J2qTJ2in n - n c{q) is the Fourier 
transform of the anomalous or Cooper channel density 
n c (x) = ?pi(x)ipi(x), where x — (x, r) comprises space 
and imaginary time labels, and the i/v are fermionic fields 
with spin label a. n c is the adjoint of n c . The remaining 
part of the action, So, which we will refer to as the refer- 
ence ensemble, describes interacting electrons in the pres- 
ence of quenched disorder, with no bare Cooper channel 



interaction. Note that even though K c = in the bare 
reference ensemble, a nonvanishing (repulsive) Cooper 
channel interaction is generated in perturbation theory 
. So thus describes a general system of disordered in- 
teracting electrons, with the only restriction being that it 
must not undergo any phase transitions in the parameter 
region we are considering, lest the separation of modes 
that is inherent in writing the action in the form of Eq. 
([?]) break down. 

We now use standard procedures [EJ to derive an OP or 
Landau-Ginzburg- Wilson (LGW) effective field theory. 
We decouple the Cooper channel interaction by means 
of a Hubbard-Stratonovich transformation. This intro- 
duces a complex valued OP field, which we denote by 
*Sf(x), that couples linearly to n c (x). The fermionic de- 
grees of freedom are then integrated out. This way we 
obtain the partition function in the form 



-F /T 



Dm 
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Here Fo is the noncritical part of the free energy, and $ 
is the LGW functional. It can be written 



$[*] = -K c / dq |tf (q)| 



In 
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where (. . .)q = J D[tp, ip] . . . e s ° denotes an average taken 
with respect to the reference action So- 

We next perform a Landau expansion or expansion of $ 
in powers of the OP. Since gauge invariance is not broken 
in the reference ensemble, only even powers of \& appear. 
The coefficients in this expansion, i.e. the vertex func- 
tions of the effective field theory, are connected correla- 
tion functions of n c in the reference ensemble. In partic- 
ular, the Gaussian vertex is determined by the pair prop- 
agator or anomalous density-density correlation function 
in the reference ensemble. Denoting the latter by C(q), 
the Gaussian term in the LGW functional reads 

$( 2 )[*] = [dq [l/\K c \ - C(q)] , (9a) 



where we have scaled \& with K c . So far we have worked 
within a particular realization of the disorder. To per- 
form the disorder average of the free enery we introduce 
replicas and use a cumulant expansion. For a particu- 
lar copy of the replicated system we then obtain, at the 
Gaussian level, Eq. ( |9a| ) with C{q) the connected prop- 
agator C(q) — {(n c (q) n c (q)) }, where {. . .} denotes the 
disorder average. C{q) is a complicated correlation func- 
tion. However, since the reference ensemble is by con- 
struction a Fermi liquid, the structure of this correlation 
function is known. Renormalization group arguments 
show that the structure of C at low frequencies and long 
wavelengths in the limit T — > is M 
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C(q) 



Z ln(rVCDq^ + |n n |)) 

h \ + {5k c /h) ln(r» /(^q 2 + |0„|)) 
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Here Oo = fc^To/?! is a frequency cutoff on the order 
of the Debye frequency (for phonon-mediated supercon- 
ductivity), and 5k c is the repulsive interaction in the 
Cooper channel that is generated within perturbation 
theory even though the bare K c vanishes. D is the dif- 
fusion coefficient of the electrons, ft. is a frequency renor- 
malization constant whose bare value is H = TzNp/2, 
and Z is the wavefunction renormalization. All of these 
parameters characterize the reference ensemble, and it is 
known that they provide a complete characterization of 
C(q) §. Using Eq. @ in Eq. (pal) we obtain 



$( 2 )[ijr] = / dq V*(q) 



t + 
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ln(O /(Dq 2 + l^nl)) 
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where t = (Sk c — Z\K C \) /h and iff has been scaled ap- 
propriately. 

Let us briefly consider this result for the Gaussian 
LGW theory. The Gaussian two-point vertex has the 
structure of Eq. (pa|), which implies Eqs. (2b). By scaling 
|q| with the correlation length £, we also obtain the be- 
havior of the latter, Eq. (fia|). Corrections to the Gaussian 

vertex are known to be of order 1/ In 2 ^f2o/(-Dq 2 + |fi n |) 

JjJ . We now need to study the higher order terms in the 
LGW functional to ascertain that they do not change 
these results, and to obtain information about the OP 
and the free energy. 

We first consider the coefficient of the l^^-term in the 
Landau expansion. It is a nonlinear anomalous density 
susceptibility in the reference ensemble which we denote 
by C^ 4 \ Due to the cumulant expansion with respect to 
the disorder average there are two different contributions 
to this coefficient, C< 4 ) = cf + cf \ where cf ] is the 
disorder average of the four-point correlation function 
for a given disorder realization, while 1S the disorder 
average of the two-point function squared. A calcula- 



tion shows that both and C^' are singular in the 
limit q — > 0. Cutting off the singularity by means of a 
wave number |p|, one finds for the leading contributions 
C[ A) ~ u 4 /|p| 4 In 4 |p| and CS[ 4) ~ V|p| 4 " d , respectively, 
with M4 and V4 finite coefficicents Q9|. The same method 
fiof shows that the most divergent contribution to the 
coefficient of the term of order \<ff\ 2n diverges like 



(4) 



C (2 



U2n 



In 2 " |p| 



(10) 



with U2n a finite coefficient. This implies that the Lan- 
dau expansion of the cutoff regularized LGW theory is 
an expansion in powers of \& /p 2 ln(l/|p|). The OP field 
theory, rather than having a simple LGW form, is thus 



strongly nonlocal. This is the technical reason for the 
unusual scaling behavior. 

The functional <!> can be analyzed by using standard 
techniques . We are looking for a FP where the func- 
tional dependence of the 2-point vertex on q and O n , Eq. 
(pc]), is not renormalized. This fixes the exponents r\ and 
z. Defining the scale dimension of the correlation length 
to be [£] = —1, Power counting shows that the coeffi- 
cients U2n (n > 2) of the non-Gaussian terms have scale 
dimensions [«2n] = (n — 1)(2 — d), and hence are irrele- 
vant operators with respect to the FP for all dimensions 
d > 2. [vi] = -2(d - 2) for 2 < d < 4, and the higher 
cumulants are even more irrelevant. The upper critical 
dimension is thus d+ = 2, and for d > 2 the critical 
behavior obtained from the Gaussian theory is exact. 

The scaling behavior of the OP, and of the free energy, 
is determined by the term of 0(\I/ 4 ) which is a dangerous 
irrelevant operator jl2| with respect to these quantities. 
For scaling purposes, the cutoff wavenumber |p| can be 
replaced by the inverse correlation length, |p| ~ 
The scaling behavior ^f> ~ p 2 ln(l/|p|) observed above 
then immediately leads to the behavior of ^f given in Eq. 
(|Ta|). One must notice, however, that the OP function ^5> 
is distinct from the physical gap function A. The latter 
determines the gap in the single-particle excitation spec- 
trum, and hence scales like the frequency or like £ 2 . This 
is expressed in Eq. ([it]). We have further corroborated 
this result by explicit calcuations. (ll| 

We next consider the critical behavior of the penetra- 
tion depth A and the conductivity a. Since we have 
shown that the mean-field/Gaussian theory yields the 
exact critical behavior at T = 0, all relations between ob- 
servables that are derived within BCS theory are valid. In 
particular, we have A ~ 1/%/A, and H C2 ~ A JI3], which 
in conjunction with A ~ £ 2 yields Eq. (^). Similarly, we 
can determine the scale dimension of the conductivity 
or resistivity. The real part of the frequency dependent 
conductivity in the superconducting phase has a singular 
contribution, which within BCS theory is given by 



Rea s (fi) = (7r 2 /2)cr n A 6(0.) 



(11) 



with a n the conductivity in the normal state. For scaling 
purposes, 5(0,) ~ I/O ~ 1/A. a n is determined entirely 
by properties of the reference ensemble, and hence it does 
not show any critical behavior and its scale dimension is 
zero. We conclude that the scale dimension of a s van- 
ishes. If we assume that the conductivity has only one 
scaling part, then the same is true for the conductivity 
or resistivity in general, and we obtain Eq. ([5]). 

We now turn to the free energy density /. Hyperscaling 
suggests that / scales like / ~ T/V ~ ^( d+2 ) j which 
leads to a homogeneity law 

f(t,T,u A ) = b-( d+ V f(t\nb,Tb 2 ,u A b 2 - d ) , (12) 

where of the irrelevant operators we have written only 
U4 explicitly. / is proportional to 114 A 4 ~ l/u4, and 
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hence the effective scale dimension of / is [/] = 4, which 
yields Eq. (Ea]). Hyperscaling is violated as usual above 
an upper critical dimension, viz. by means of the quartic 
coefficient being dangerously irrelevant with respect to 
the free energy. By differentiating twice with respect to 
T one obtains the specific heat coefficient, Eq. (f4~E|). 

The quartic term whose coefficient is yields cor- 
rections to scaling that represent fluctuations in the po- 
sition of the critical point: By making K c a random 
variable, and integrating out that randomness, one ob- 
tains a term of that structure. Repeating the arguments 
of Ref . H , we find that the relative fluctuations of the 
position of the critical point decay anomalously slowly, 
At/t ~ £-( d ~ 2 ). (5-correlated elastic scatterers thus 
induce long-ranged correlations of the t- fluctuations. 
Translating that into the corresponding fluctuations of 
T c via T c = T exp(-l/|i|) we obtain Eq. §). 

Let us now discuss the physics behind the unusual crit- 
ical behavior that we have found. The reason for the non- 
local LGW theory and the ill-behaved Landau expansion 
is that in deriving the OP field theory one integrates out 
the fermions. This procedure leads to a local OP theory 
only if the OP fluctuations or critical modes are the only 
soft modes in the system. This is not necessarily the case, 
there may be soft modes that are distinct from the OP 
fluctuations, but couple to the latter. These modes may 
be due to conservation laws, broken symmetries, or long- 
range forces. This caveat holds for all phase transitions, 
including thermal ones. However, at T — the number 
of soft modes is larger than at finite T, and therefore a 
nonlocal OP theory is more likely. Here the additional 
soft modes are particle-hole excitations in the Cooper 
channel that are soft only at T = 0. They are integrated 
out in deriving the OP functional, and they lead to both 
the logarithmic structure of the Gaussian vertex, and to 
the divergence of the higher vertices. While this leads to 
a breakdown of the LGW theory, it also is the reason why 
the critical behavior can be determined exactly despite 
the ill-behaved nature of the OP field theory: The addi- 
tional modes lead to an effective long-range interaction 
between the OP fluctuations (one that falls off only log- 
arithmically at large distances and times). For thermal 
phase transitions it is well known that long-range inter- 
actions stabilize mean-field theory ]l5j , and the same is 
true here. Notice that in the present case the long-range 
interaction is self-generated by the system. 

The same basic physics as discussed above, namely soft 
modes other than the OP fluctuations, also makes the 
quantum critical behavior of itinerant ferromagnets non- 
trivial, and yet exactly soluble p6| . There are, however, 
important differences between the two phase transitions. 
In the magnetic case, the additional modes are weaker, 
which leads to a Gaussian FP with standard power-law 
critical behavior. Furthermore, several competing dy- 
namical or temperature scales lead to complicated dimen- 



sionality dependencies of critical exponents even though 
the FP is Gaussian. In the present case the additional 
modes are so strong that they completely dominate the 
physics. As a result, there is only one temperature scale, 
with a scale dimension [T] = 2, and the quantum critical 
behavior of disordered bulk superconductors is BCS-likc. 
However, strong disorder fluctuations lead to a broad 
transition region as d — ► 2. 
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